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ABSTRACT 

A theory of pack aluminization haj been formulated by combining 
gaseous and solid-state diffusion rates. This theory relates the 
surface composition of the coating and therefore, in principle, the 
phase morphology and the growth rate of the coating, to pack operating 
parameters such as pack A1 density, type of activator, temperature and 
ot!iers. Experimental data on the aluminization of unalloyed Ni In 
pure A1 packs obtained by us to date are in good agreement with the 
predictions of the theory. 
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ABSTRACT 

A theory of pack aluminization ha: been formulated by combining 
gaseous and solid-state diffusion rates. This theory relates the 
surface composition of the coating and therefore, in principle, the 
phase morphology and the growth rate of the coating, to pack operating 
parameters such as pack A1 density, type of activator, temperature and 
others. Experimental data on the aluminization of unalloyed Ni in 
pure A1 packs obtained by us to date are in good agreement with the 
predictions of the theory. 
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ON THE KINETICS OF PACK ALUMINIZATION 


Introduction 

The pack aluminization process Is extensively used to form pro- 
tective coatings on superalloy components of aircraft gas turbines. 

In this process the metal to be coated Is placed In a powder mixture 
consisting essentially of Al* a halide activator such as NaF* and an 
Inert filler material such as AI 2 O 3 . Other metals such as N1, Cr, 

SI, T1 are sometimes added (1, 2) to control the A1 activity or other- 
wise Influence the properties of the coating. During aluminization, 

A1 Is transferred from the pack to the specimen surface In the form of 
aluminum halides. These compounds decompose at the surface to release 
A1 which diffuses Into the metal, creating the Al-rich coating. The 
aluminization process thus Involves three steps: 1) gaseous diffusion, 

2) surface reaction and 3) solid diffusion. Of these three it Is 
assumed that the surface reaction step Is rapid and the kinetics of the 
process Is governed by a combination of gaseous diffusion in the pack 
and solid state diffusion in the coating. The gaseous diffusion step 
was recently analyzed by Levine and Caves (3). In the following a quant- 
itative description of the kinetics of aluminization Is given which takes 
diffusion in the solid as well as the gas phase Into account. 

Kinetics of Pack Aluminization 

A. Gaseous Diffusion : Levine and Caves' (3) model which applies to 
packs containing unalloyed Al, assumes that aluminum is transferred from 
the pack to the coating surface by the diffusion of aluminum halides 
through an Al-free zone in the pack adjacent to the specimen surface. 



In the body of the pack, the halide vapor pressures are taken as those 
in equlllbrlinn with the bulk pack composition. The halide vapor pressures 
at the specimen surface, which is assumed to be at a constant A1 activity, 
are determined by mass-balance and thermodynamic equillbriian conditions 
at this surface. Chemical reactions are assumed to occur rapidly enough 
at the pack-coating interface to allow the establishment of a local 
equilibrium. 

As A1 is transformed from the pack to the specimen surface it is 
assumed that the width of the depleted zone and consequently the dif- 
fusion distance for the halide vapor Increases. This leads to a para- 

2 

bolic relationship between the weight of A1 transferred Wg, in gms./cm 
and time, t, in seconds. 


y2 , 2£eM Al__ ^ a ^ 
g i g 


( 1 ) 


in which e and i are correction factors for pack porosity and pore length, 
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M is the atomic weight of A1 , p is the pack A1 concentration in gms/cm , 

2 

d is the diffusion distance in cm, A is the total area in cm , and 
is the overall rate of diffusion of A1 through the gas phase in moles/sec. 
The instantaneous flow of A1 through the gas phase is given by 


Uai£ ■ 1 £ D, (P,-P;) (2) 

ITT 1 ’ ' ' 


where is the diffusivity, and are the partial pressures of the 


i— Al-bearing species in the pack and at the coating surface, respectively, 
T is the temperature and R the gas constant. Hence the parabolic rate 
constant. Kg, is given by 



( 3 ) ' 


Application of the above equations Involves calculation of the 
equilibrium pressures and P^. The method of calculation depends 
upon pack operating conditions and type of activator used. 

The pack aluminization process Is generally carried out at atmo- 
spheric pressure by placing a semi -open retort with the pack and the 
parts to be coated In a flowing H 2 or Inert atmosphere. The activators 
In common use can be divided Into two classes. The sodium halides 
1 .e.,Na(F, G, 6 r, I) as well as AIF^ remain In the pack as condensed 
phases at the usual activator ratios (l-io'^/o) and pack cementation 
temperatures (800-1100°C). If NH^F or NH^F.HF are used as activators, 
crystalline AIF 3 forms In the pack by reaction with aluminum. With 
the highly volatile activators A 1 (Cl, Br, I)j and NH^{C1, Br,I), no 
condensed phase Is formed In the pack. The composition of the gas- 
mixture that forms In the pack will be unstable because of continuous 
losses. While the function of both kinds of activators Is qualitatively 
Identical, only the former class can be quantitatively analyzed. The 
composition of the gas phase can be calculated given the temperature. 
alumlniOT activity and the total pressure. For NaX, the partial pressure 
ip^) of the significant species In the bulk pack will be given by the 
equilibrium constants of the following chemical reactions, subject to 
the last two conditions: 


NaX [i) + A1_ ■ AlX(g) + Na(g) 

(4) 

2 NaX (A) + ^ - AlX 2 (g) + 2Na(g) 

(5) 

3 NaX U) + ^ • AlX 3 (g) + 3Na(g) 

( 6 ) 

NaX il) - NaX (g) 

(7) 


NaX(l) + h H 2 (g) • HX(g) ♦ N<(g) (8) 

2 A1X3(g) • Al2X5(g) (9) 

'’N.(g)* 3 ’’aIXj,,)*® \Xg(9f 2 '’a1X2(,) 

* ’’AlX(g)* ’’HX(g) (10) 

z * 1 atm. (11) 


For AlFj or ammonium fluoride activators, the partial pressure of AlFj(g) 
equals the vapor pressure of AlFj(g) and Is Independent of the aluminum 
activity. The other partial pressures can be calculated from: 

AlF 3 (g) + 2^-3 AlF(g) (12) 

2AlF3(g) 3 AlF 2 (g) (13) 

AlF 3 (g) + |H 2 (g) » 3 HF(g) + M (14) 

’’AlF 3 (g) * '’AlF 2 (g) '*HF(g) ^H 2 (g) “ ^ 

Similar calculations for at the specimen surface must take Into 

account certain kinetic considerations. The pack gases not only trans- 
port Al, but also halogen, hydrogen and sodium gas. For the NaX activators, 
Na(g) and halogen must be transported In equimolar proportions, so that 
they react and condense back Into NaX(i) near the specimen surface. 
Therefore, eq. (10) Is modified to: 


^Na^'^Na ‘ ’’fia^ “ ^ *^ 3 ^^A 1 X 3 " '’^1X3^ * ^^1X3^ * ^I^^AIX" ^AIX^ 


°HX^^HX ■ ^HX^ ® °6^’’a 12X5 ‘ ^Al2Xg^ 


HX^ 


(16) 


For AIF 3 , since Its vapor pressure must be constant throughout the pack, 
a part of the Al that Is carried by AIF and AIF 2 to the specimen surface 
Is lost due to formation of AlF 3 (s), for example: 3AlF(g) « AlF 3 (s) + 2 Al. 
This Implies that the actual Al transfer equation should be: 


W. 


[{2/3).D^.(P^,P 


- Paif) 


+ (1/3).D,.(P 


AlF ■ ^AIF 


)].t 


( 17 ) 


In principle, a H 2 balance equation should be Incorporated, but has 
been found to be practically unimportant. The A1 deposition rate con- 
stant, Kg, has been computed at SOO^C and 1093°C for NaC1, NaF, Nal 
and AIF^ as a function of specimen surface composition. The A1 activity 
as a function of composition In the N1-A1 system (4) and thermochemical 
data (5) are available In the literature. The gas diffusivltles and 

the value of t « 4 were taken from Levine and Caves' paper (3). The 

3 

pack compositions we selected corresponded to 0 > 0.04 gm/cm and 
porosity, e » 0.70. 

B. Solid Diffusion ; If the surface of the specimen remains at constant 
composition, assuming that diffusion takes place Into an Infinitely ex- 
tended medium, the rate of diffusion In the solid will also obey a par- 
abolic law expressed by the equation 

W5 = Kjt (18) 

where K^, the parabolic rate constant for diffusion In the solid, will 
depend upon the value of the surface composition. It is evident that 1n 
order for the con^ositlon at the surface of the coating to be time in- 
variant, the K values must be identical, I.e., the surface composition 

must assume that value for which ® K„. This unknown value of surface 

s g 

composition can be predicted if K^, as well as Kg, can be calculated. 

For complex superalloys, the necessary diffusivity and thermodynamic 
data are not yet available to carry out the calculation of K^. Sufficient 
data are available for the binary N1-A1 system, however, to permit calcu- 
lation of the Kg values for the aluminization of unalloyed N1. 


RKPROBUCmrUTY OF THE 

' ' : IS POOR 


Hie aluminide coating on nickel will ctmsist of single*‘phase 
intermetallic layers of progressively decreasing A1 content, such as 
Ni^Al^, NiAl, Ni^ and the Ni(Al) terminal phase. The id^tity of 
the first layer is determined by the surface con^wsition. The position 
of the interface between any two layers i and j is given by: 

The free surface of the coating also moves outward parabolically with 
time, because A1 is gained but Ni not lost in pack aluminization (6) . 
The thickness of the layer j is: 


tj “ CPij - Pjjj) » qj/t C20) 

The growth rates of the various layers, as well as of the coating as a 
whole are known if the constants p^^ can be calculated. These dep«id 
on the interdiffusion coefficients and the composition parameters of 
the various phases. 

The first step in calculating p^^ is to set 15) material balance equations 
at the interfaces. Thus, the rate of movement of the ij interface is 
given by: 


■ -DiOCi/3x)y * Dj{3Cj/5x)j. 

Cj^j and Cj^ are the equilibrium compositions, in atom fraction, of the 
phases i and j, respectively, at the interface ij. For a surface con^sosition 
Cg, the material balance equation for the free coating surface i.e., the 
vapor-solid interface, s, is: 


(1 - Cg) (dXg/dt) - (3c^/3x)g 


(21b) 


If the interdiffiision coefficients do not vary with ccmiposition, then 
analytical solutions are known (7,8) for the composition profile in each 
phase, for instance the solution for the phase j is: 



( 22 ) 

With the help of these and eq. (19), all tte derivatives in eq. (21a,b) can 
be replaced, and a system of n equations in n unloiowns , p^^^ , idwre n is 
tl» total nunber of phases In the coating, can be generated and solved. 

This analytical scheme suffices in many practical cases of binary dif> 
fusion. In the Ni-Al system, however, it has been found that idiile 
Ni 2 Alj, Ni^ and the Ni(Al) texminal solution have diffusivities 
practically invariant with conposition, the interdiffusion coefficient 
in KiAl varies strongly with cm^sition (9) . It has a minimun at 
around the equiatomic con^sition, and oi either side of this con^sition, 

D varies exponentially: 

'VfiAl ■ ^ C®'’ (23) 

For such a case, an analytical solution is not known. Therefoi^, the 

con^sition profile in NiAl has to be ccmpited numerically. This is 

done by first transforming the diffusion equation in NiAl: 

(3c/3t) » 3/3X P (3c/3x)) (24) 

to an ordinary differential equation: 

y(dVdz^) ♦ (P 12/2 Dq) 2 Cay/dz) • 0 (25) 

by substitution of the new variables: 

y ■ C3q) (8c) and z ■ (x/p.,/E) [p,,« rate const, for the 

NiAl/NijAl interface 

The transfoimed eq. (25) can be solved by a numerical procedure described 
in detail in Ref (8). In Fig. 2, we show a number of concentration pro- 
files across NiAl-NljAl-Ni(Al) terminal solutirai type of coatings, num- 
erically confuted for 20 hours of pack cementation at 1100*C. The 
numerical method becomes important only when NiAl is the surface layer 


c.(x,t) ■ c.. 


erf(x/2vO) - erf (Pjj/2»^ 




of tht coating. If th« surface composition Is such that NI 2 AI 2 Is the 
first layer* then we have found that all the other layers develop relatively 
small thicknesses. In these cases* a constant value of the Interdiffusion 
coefficient for N1A1, obtained by averaging the D over Its homogeneity 
range* can be used In the analytical scheme to give accurate enough 
results for NlgAlj-NIAl-NIjAl-NlCAl) type of coatings. 

Thus* by the appropriate application of the analytical or numerical 
method the growth*rate constants p|j or and the composition profile 
across the coating C 2 n be calculated* for any surface composition. The 
A1 Intake of ^e nickel specimen* W^* and the rate constant are then 
obtained by Integrating the profile. 

In Fig. 1 (a) and (b) the versus surface composition curves at 
1093^0 and SOO^C have been superimposed on the curves to locate their 
Intersections. 

Comparison with Experiment 

Data have been obtained for the aluminization of unalloyed N1 in 
packs activated with AIF^* NaF* NaCI and Nal. The packs contained 

4 '^/o AT * 4 ^/o activator* balance AlgOj* and were contained In alumina- 
lined Iron retorts with a slide-fitting cover. The retorts were heated 
In an atmosphere of After removal* the coated specimens were cleaned 
by brushing and In an ultrasonic bath. Surface compositions were determined 
by microprobe analysis of the as-coated surfaces. 

The square of the weight gain and surface composition vs. time curves 
for a nimiber of specimens are shown In Fig. 3. It was observed that 
surface melting had occurred on specimens aluminized at 1093°C In packs 


T> ^ > 





activated with AIF^ and NaF> and curves for these specimens are not 
given. Also low and erratic weight gains were obtained with Nal ac- 
tivated packs* and these results are not shown. For the t'emalnder of 
the activators and ten^eratures studied It was observed that the surface 
compositions of the coated specimens were relatively Independent of time. 
Furthermore the square of weight gain vs. time curves were nearly 
straight line obeying parabolic relationship. These results Indicate 
that within a short period the system reaches steady-state at the pack 
coating Interface. The surface con^osltlons obtained under various 
ccidltlons are cons^ared with those predicted theoretically In Table 1. 

It will be noted* the activators arc properly ordered by the theory. 
Furthermore* there Is good correspondence between the predicted and 
observed values of surface composition. This is felt to be a strong 
confirmation of the basic validity of the proposed model. 

Summary and Conclusions 

By confining the gaseous diffusion model of Levine and Caves with 
a calculation of diffusion rates In the solid a theory of pack alinnlnl- 
zatlon has been formulated. This theory relates the surface conqsositlon 
of the coating (and therefore, In principle* the coating formation rate) 
to pack operating parameters such as pack A1 density* type of activator* 
temperature and others. Experimental data on the alimilnlzatlon of un- 
alloyed N1 In pure A1 packs obtained by us to date are In good agreement 
with the predictions of the theory. 


Acknowledgement: This work was supported In part by the NASA Lewis' 
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Table 1. Comparison of theoretically predicted and observed results 


Pack Conditions 

1 

Phase; Surface con^)cs1t1on (a/o Al) 1 

1 

Activator 

Coating Ten^rature 

Observed 

Predicted 

AIFj 

1093°C 

Liquid, >62.5 

Liquid, >62.5 

A 1 F 3 


NI 2 AI 3 , 59.0 

NI 2 AI 3 , 59.10 

NaF 

1000°C 

NIjAlj 

NIjAlj ■ 

NaCI 

1093®C 

N1A1, 53.2 

N1A1, 54.0 
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